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From an optical spectroscopic survey of 3CR radio galaxies with z < 0.3, we discovered a new spectroscopic class of powerful radio- 
loud AGN. The defining characteristics of these galaxies are that compared with radio galaxies of similar radio luminosity they have: 
a [O III]/U/3 ratio of ~ 0.5, indicative of an extremely low level of gas excitation; a large deficit of [O III] emission and radio core 
power. We interpret these objects as relic AGN, i.e. sources that experienced a large drop in their level of nuclear activity, causing a 
decrease in their nuclear and line luminosity. This class opens a novel approach to investigating lifetimes and duty cycles of AGN. 
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1. Introduction 

Optical spectroscopy has played a major role in enhancing our 
unde rstanding of active g alactic nuclei ( AGN~) . iHeckmanl (1 1 980h 
and lBaldwin et al.l d!981l) demonstrated that optical lines can be 
used as tools to classify in general emission-line objects, particu- 
larly AGN. Diagnostic diagrams comparing emission line ratios 
can distin guish H II regions from gas clouds ionized by nuclear 
activity dVeilleux & Osterbrocklll987D . Moreover, AGN can be 
separated into Seyferts and Low Ioniz ation Nuclear Emission- 
line Regions (LINERs. lHeclmianll980h_ since they form separate 
branches in the diagnostic diagrams dKewlev et al. 2006). 

We performed an optical spectroscopic survey of the 113 
radio gal axies (RG) belongi ng to the 3CR sample and with 
z < 0.3 (Spi nrad et alJ 1985h. using the Telescopio Nazionale 
Galileo (Buttiglio ne et al.l 120091 1201 Obh . Most RGs belong to 
two main spectroscopic classes, those of high and low excita- 
tio n galaxies (HEG a nd LEG respectively, originally introduced 
the analogous to Seyf ert and LINER for 



bv lLaingetal]|1994 



radio-loud AGN). In lButtiglione et ali d2010al) . we also reported 
the discovery of a new class, characterized by an extremely low 
level of gas excitation whose properties are discussed in this 
Letter. 

We adopt the following cosmological parameters: H — 71 
km s _1 Mpc" 1 , Q A = 0.73, and Q m = 0.27. 



Send offprint requests to: A. Capetti 

* Based on observations made with the Italian Telescopio Nazionale 
Galileo operated on the island of La Palma by the Centro Galileo Galilei 
of INAF (Istituto Nazionale di Astrofisica) at the Spanish Observatorio 
del Roque del los Muchachos of the Instituto de Astrofisica de Canarias. 



2. A new spectroscopic class: 'Relic' radio galaxies 

Fig.Q]shows the location of the 3CR sources in the emission-line 
diagnostic diagrams. They all fall in the area populated by AGN, 
with only one exception represented by an object with emission 
line ratios typical of star-forming galaxies. Their [OIITJ/H/? ra- 
tios are in the range 1.3 - 16. 

Three galaxies (namely 3C 028, 3C 314.1, and 3C 348, 
shown in Fig. [2] as exampl^B) s tand out for a lower ratio, 
[O III]/Hy6~ 0.5. For this reason, in lButtiglione et al.1 (1201 Oal) we 
refer to this class as 'extremely low excitation galaxies' (ELEG). 
They are well separated from the rest of the 3CR sample and 
fall in a region scarcely p opulated by Sloan Digital Sky Survey 
(SDSS, iKewlev et al.ll2006l) emission-line sources. Considering 
first their [S II]/Ha' and [O l]/Ha ratios, they are located in 
the region of the diagnostic diagram occupied by AGN. In the 
[OIII]/Hy6 versus [N II]/Ha' plane, they straddle the boundary 
separating AGN from composite galaxies, which are possibly 
objects with a contribution from young stars and AGN to the 
emission lines. 

The ELEG are also distinctive because of their radio proper- 
ties (see Table. [TJ: 3C 028 and 3C 314.1 do not appear to con- 
tain a radio-core at the level of 0.2 and 1.0 mJy respectively (see 
iGiovannini et al.|[l988l and refere nces therein), while a 10 mJy 
radio core is detected in 3C 348 (Mor ganti etal1ll993h . The ra- 
tio of core (at 5 GHz) to total (at 178 MHz) radio emission for 
these sources is in the range F core /F tot < 10 5 - 10~ 4 , compared 



1 Since the [S II] doublet is not covered by our spectrum, 3C 348 
is not represented in the central panel of Fig. [TJ where, for the same 
reason, the star-forming galaxy discussed above is also missing. 
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Fig. 1. Contours: density of SDSS emission-line galaxies in the optical diagnostic planes (adapted from iKewlev et all 12006). 
Galaxies below the curved solid lines are star-forming galaxies. In the left diagram, the dashed curve marks the transition from 
composite galaxies to AGN; in the middle and right panels, the oblique line separates LINERs from Seyferts. The circles are the 
3CR sources. These diagrams reveal the existence of a class of extremely low excitation galaxies, i.e. the three filled red circles 
corresponding to the objects with the lowest [OIII]/Hy6 ratios. 



Table 1. Properties of the relic radio galaxies 



Name 


z 


L[OIII] 


Lcore 


Ll78 






erg s _1 


erg s" 1 Hz 1 


erg s Hz 


3C 028 


0.195 


40.96 


<29.33 


34.24 


3C 314.1 


0.120 


39.69 


<29.56 


33.59 


3C 348 


0.154 


40.41 


30.80 


35.35 



to an average ratio of ~ 1 .6 x 10~ 3 for RGs of similar luminosity. 
They also have a [OIII] line luminosity that is 10 to 1000 times 
lower than RG of similar radio power (see Fig.00 

On the basis of these findings, we interpret these objects as 
relic AGN in which the nuclear activity is currently switched- 
off (or substantially weaker than its long-term average level). 
The quantities considered (radio-core, emission-line luminosity 
and ratios, and total radio emission) respond to changes in the 
AGN activity level on different timescales. This is because they 
are produced on scales of sub-parsec for the radio-core, to kpc 
for the line emission produced in the narrow line region (NLR), 
to hundreds of kpc for the extended radio-emission. A drop in 
the activity level is rapidly followed by a decrease in both the 
radio core flux and the flux of ionizing photons. This in turn 
causes the cooling of the NLR, leading to a relative low level of 
[OIirj/15007 emission and a decrease in the [OIII]/H/3 ratio, as 
described in more detail in the next section. In contrast, the lu- 
minosity of the extended radio emission remains essentially un- 



2 As noted above, two ELEG lie in Fig.[T]among the composite galax- 
ies, generally thought to be objects with a mixed contribution from 
young stars and AGN to the emission lines. However, if this were the 
case for the ELEG, they should exhibit a line excess with respect to RG 
of similar radio luminosity, in contrast to what is observed. 



changed over very long timescales, accounting for the observed 
properties of ELEG. 



3. Spectroscopic time evolution of relic RG 

Binette & Rob inson] ( 1 19871) investigated the time evolution of the 
properties of a photoionized cloud following a sharp decrease 
in the strength of the ionizing photon field. They considered a 
plane-parallel cloud, in the low density limit, and an instanta- 
neous drop in the intensity of the nuclear emission by a factor of 
100. Their analysis found a very rapid quenching of the [OIII] 
line. This is due to the high efficiency of the charge exchange 
reaction, + 2 + H° — > + + H + , which proceeds rapidly as 
soon as the gas begins to recombine and the fraction of neutral 
hydrogen increases. They estimated that the decay time of the 
[O III] (the time after which the line intensity has decreased by 
a factor 1/e) is f^([0 III]) ~ 2000 rC e years, where n e is the elec- 
tronic density. Since the decay time of the Hfi line is substantially 
longer, t d (H/3) ~ 130000 n; 1 years, the [O HTJ/H0 ratio decreases 
with time. Starting from an initial value of [OIII]/HjS~15, after 
~ 4000 n^ 1 years the ratio has decreased to ~ 0.5 (the value ob- 
served in ELEG). 

All relevant timescales depend on density. The density trac- 
ers corresponding to the brightest emission lines are the [O 
II]/L13726,3729 and [S II]/U67 16,6731 doublets. The individ- 
ual measurement of the [O II] lines is irreparably compromised 
by the intrinsic line widths in the NLR, leaving one only with the 
possibility of exploiting the [S II] doublet alone. Unfortunately, 
from our spectra we cannot obtain an accurate estimate of this 
parameter. We have no measurements of the [S II] for 3C 348 
(not covered by the TNG data, see the right panel in Fig. 0, 
while in our spectrum of 3C 028 the [S II]/16731 line falls in a 
CCD defect. Only for 3C 314.1 are we able to measure a ratio 
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Fig. 2. TNG spectrum of 3C348, an example of a relic RG. The left and right panels are centered on the H/3 and Ha spectral regions, 
respectively. The original spectra are indicated by the upper solid lines with superimposed we show the best-fit stellar populations 
model used to subtract the starlight; the lower solid lines are the residual spectra. The flux is given in 1CT 18 erg cirT 2 s _1 A , while 
the wavelength is in A, in the rest frame in the lower axis and in the observed frame in the axis above. The three main telluric 
absorption bands are indicated as circles with a cross inside. 



[S H]/16716/[S IIM6731 = 1.6 ± 0.3, whic h is indicative onl y of 
a low density (n e < 3 x 10 2 cirT 3 ) regime dOsterbrocklll989l) . 

An additional factor to consider in the temporal evolution is 
the light travel time from the nucleus to the NLR clouds. The 
light travel time from the nucleus to the edge of the 2" slit is 
~ 10 4 years. This must be considered as an upper limit to the re- 
sponse of the portion of the NLR probed by our observations to 
a change in the nuclear flux since emission lines of R G are usu- 
ally h ighly concentrated around the nucleus (e.g. lTremblav et all 
2009). Nonetheless, considering the uncertainties related to the 
actual distribution and density of ionized gas, it can be envisaged 
that the light travel effects might play a key role in the evolution 
of the observable properties of the NLR. 



4. Discussion and conclusions 

The energy carried by the jets of radio-loud AGN has a profound 
impact on the evolution of the ir hosts and the energy balance 
of the intracluster medium (e.g. Fabian et al. 2003: Croto n et al.l 
2006). To ascertain the effects of AGN feedback, it is necessary 
to measure not only the jet power, but also the timescales and 
frequency of the energy release, i.e. measure AGN lifetimes and 
duty cycles. Several approaches, base d purely on radio proper- 
ties, have been used previously (e.g . lAlexander & Leahvll 19871: 
Shab ala et ai]l2008l: iBird et al.1120081) to investigate these issues 
These analyses ha ve led to evidence dGiovannini et al.l fl988t 
Parma et al] 120071) of relic (fossil) RG based on the low core 
dominance or very steep spectral indices (a manifestation of the 
energy losses of the radio-emitting relativistic electrons) of some 
systems. Using radio observations it is possible to derive model- 
dependent estimates of the RG lifetimes or, in case of fossil RG, 
the time elapsed since the 'death' of the active nucleus, tq. 

Relic RGs, distinguished by their optical spectroscopic prop- 
erties, open a new and complementary path to exploring the life- 
time and duty cycles of radio galaxies. The fraction of relics in 
the RG population is given by the ratio of the duration of the 
relic phase to the RG lifetime. In principle, both parameters can 
be derived once all RG in a complete sample are spectroscopi- 
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Fig. 3. [OILTJ vs. total radio luminosity at 178 MHz (in units of 
ergs ' and ergs 1 Hz 1 respectively) for the 3CR sources with 
z<0.3 (limiting for clarity to its bright luminosity end). HEG are 
marked with circles, LEG with squares. The green asterisks are 
ELEG, i.e. the relic RG. The triangles are objects that could not 
be classified spectroscopically. The dotted lines are the best-fit 
linear relation between L[oni] and Ly]% obtained for HEG and 
LEG separately. 



cally classified, a situation that is (almost) met by our 3CR spec- 
troscopic survey. 

Considering the range of extended luminosity where relic 
RG are observed, Log L\j% > 33 ergs 1 Hz , there are 83 ob- 
jects in the 3CR with z < 0.3, three of which are relics, one has 
a star-forming-like spectrum, 46 are HEG, and 20 are LEG. This 
leads to the fraction of relics being 3/83. This fraction can in- 
crease substantially depending on the nature of the 13 objects in 
the sample that could not be classified spectroscopically. In par- 
ticular, as explained in more detail below, some of them can be 
considered as 'candidate relics'. A second, less important, un- 
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certainty in the relic fraction would also exist if HEG only can 
evolve into relic systems. 

The duration of the relic phase could in principle be esti- 
mated using the spectroscopic evolution, by comparing the tem- 
poral changes in line ratios with those predicted by photoioniza- 
tion models. Furthermore, the intensity of the various lines can 
be compared with those of RG of similar extended luminosity, 
i.e. with the relic progenitors. This approach requires a direct 
measurement of the NLR density, the parameter that drives its 
temporal evolution. Unfortunately, the density in the relic galax- 
ies is essentially unconstrained by our observations as we are un- 
able to measure the [S II] doublet ratio in two of the ELEG, and 
it is poorly constrained for 3C 314.1. However, future measure- 
ments should be easily attainable with dedicated observations as 
the intensity of the [S II] doublet is similar to that of the Ha line, 
whose flux was accurately measured in our spectra. 

Although the extended radio emission is, in this context, the 
structure responding on the longer timescale to the changes in 
the nuclear properties, the radio morphology of the three relics 
can provide useful insights into their evolution. All of them 
can be classified as FR II following the original definition of 
iFanaroff & Rilevl d 1974b . but they differ when examined in more 
detail: 

- 3C 028 has well-defined twin jets linking the host to a 
double-lob ed structure, with w ell-defined hot spots, but lacks a 
radio core ( Feretti et al . 1984). There is no apparent brightness 
discontinuity in the jets on either side of the host, an indication 
that the drop in nuclear activity occurred recently. 

- 3C 3 14. 1 has a structure known as a "fat double" as it lacks 
jets and hot spots dLeahv & Perlevl[l99ll) . This is an indication 
that the switching-off of the jets of this galaxy occurred at least 
~ 7 x 10 5 years ago, based on the light travel time to the edge of 
the radio source ~ 250 kpc from the nucleus. 

- 3C 348 is considered to be a "born again" RG based on 
its peculiar radio morphology dGizani & Leahvl l2003): two jets 
emerge from the core and propagate within a relaxed double- 
lobed, steep-spectrum structure. On the basis of the size of the 
inner radio structure embedded in the FR II fossil (~ 200 kpc), 
the relic phase started at least 3 x 10 years ago. Despite the 
presence of a restarted AGN, the relic line emission appears to 
be still the dominant component observed in this galaxy; this 
implies that the current phase must be associated with a signifi- 
cantly lower level of nuclear luminosity than the earlier phase of 
activity. 

Evidence of the different evolutionary stages in the three 
sources is supported by comparing their [O III] luminosity with 
respect to HEG of similar radio luminosity (see Fig. |3j: for 
3C 028 this is lower by only a factor of < 10, while 3C 314.1 
and 3C 348 exhibit much larger deficits (a factor of 10 2 - 10 3 ) . 

An independent estimate of the relic phase duration, Td, can 
be obtained by directly spatially mapping the change in state of 
the AGN using deeper o ptical spectroscopy of the off-nuclear 
emission line regions (e.g. Robins on et al.| [l987). Owing to light 
travel effects, information about the drop in the nuclear emission 
might not have reached these outer regions. In a sort of reversed 
light-echo effect, the regions located at radii r > ctu would 
still show the original high excitation state, characterized by a 
ratio [OIII]/H/J ~ 20 times higher than in the nuclear regions. 
In practice, geometrical effects affect the estimate of td, since 
the isochrones are paraboloids centered on the nucleus and the 
derived timescale depends on the orientation of the gas cloud 
with respect to the line of sight. This degeneracy can be broken 
when emission lines are detected on both sides of the nucleus. 



Finally, we consider in more detail the 13 spectroscopically 
unclassified sources (usually because the [O III] and/or the H/3 
lines are too faint to be measured). Not all of them should be 
considered as potential additional relics since the uncertain clas- 
sification is caused by a variety of factors^. However, relative 
to RG of similar total radio luminosity nine unclassified sources 
also have a relatively low [OIII] luminosity (see Fig. [3} and we 
therefore consider them to be plausible 'candidate relics'. The 
importance of these objects is two-fold: first of all, only after a 
proper spectroscopic validation will it be possible to derive the 
ratio of the number of relic to active galaxies, and consequently 
the relative duration of the two phases; secondly, it is possible 
that these sources might be in a different evolutionary stage than 
the three relic galaxies discussed here. A clearer characterization 
of the class of relic RG would require a study of these objects 
in greater detail from both the radio and spectroscopic point of 
view. 

Relics AGN are also likely to exist among radio-quiet ob- 
jects; for example, NGC 5252, with its high-excitation extended 
NLR surrounding a nuclear region with a LINER spectrum, i s 
a likely relic QSO (iGoncalves et al.lll998t ICapetti et al.l l2005). 
In the region of the spectroscopic diagrams typical of the 3CR 
relics there is indeed a substantial number of SDSS emission- 
line galaxies (see Fig. 1). However, the lack of radio diagnostics 
makes it identify genuine relics among the broad distribution of 
radio-quiet AGN in the diagnostic planes. Furthermore, relics 
are likely to be rare in a flux-limited sample of emission line 
galaxies, particularly when their definition requires the detection 
of the short-lived [O III] line. As an alternative, one could rely 
on non-classical diagnostic diagrams, including lines originating 
from high excitation gas unaffected by the charge exchange reac- 
tion (e.g. Hell/14686) that causes the extremely short-decay time 
for the [O III] line. The temporal evolution within these same 
diagnostic diagrams can provide us with an estimate of the relic 
age. However, the likelihood of observing these effects will in 
general depend on both the gas density and (through light travel 
times) the size of the emitting region. In principle, extrapolating 
backward in time it will be possible to associate these sources 
with their active counterparts at the correct level of luminosity 
and thus estimate their lifetimes. 
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